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ABSTRACT 



Aims. We attempt to understand the presence of gas phase CO below its sublimation temperature in circumstellar disks. We study two 
promising mechanisms to explain this phenomenon: turbulent mixing and photodesorption. 

Methods. We compute the chemical evolution of circumstellar disks including grain surface reactions with and without turbulent 
mixing and CO photodesorption. 

Results. We show that photodesorption significantly enhances the gas phase CO abundance, by extracting CO from the grains when 
the visual extinction remains below about 5 magnitudes. However, the resulting dependence of column density on radial distance 
is inconsistent with observations so far. We propose that this inconsistency could be the result of grain growth. On the other hand, 
the influence of turbulent mixing is not found to be straightforward. The efficiency of turbulent mixing depends upon a variety of 
parameters, including the disk structure. For the set of parameters we chose, turbulent mixing is not found to have any significant 
inffuence on the CO column density. 
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1. Introduction 

In circumstellar disks, because of the high densities, carbon 
monoxide is expected to stick onto grains at temperatures below 
about 17 K. However, millimeter interferometric observations of 
' several circumstellar disks around T Tauri stars (e.g. DM Tau, 
', GM Aur, and LkCal5) have revealed that a large amoun t of CO 
■ remains gaseous at tem peratures as low as 10 K (Dai'to is~et alj 
. ^ ; 120031; |Pi6tu et alj|2 007). These authors show that the bulk of CO 
, is at temperatures below 17 K between 100 and 800 AU from 
5_j the central star Since these disks are much older than the stick- 
Ci ing timescale of CO onto grains (~ 10^ yr), some mechanism 
must prevent the complete depletion of CO. 

lAikawa & Nomural (l2006h proposed that the presence of 
cold CO in disks can be a conseque nce of ve rtical turbulent 
mix ing (See also Fig. 4 of I Willacy e t al. 2006). ISemenov et alj 
( 1200 6) studied in detail the influence of both vertical and radial 
mixing on the abundance of cold CO. They found that vertical 
mixing cannot account for the observed abundance of cold CO, 
while the situation improves slightly when considering the full 
2D mixing (vertical and radial). However, Aikawa (2007) com- 
puted analytically the influence of vertical mixing on the abun- 
dance of cold CO and concluded that vertical mixing was suf- 
ficient to achieve an abundanc e of CO consistent with that ob- 
served by Dartois et akj (l2003h . The apparent discrepancy with 
[Semenov et al. (2006) can be understood as a consequence of 
the different effici encies of grain surface reactions in their stud- 
ies (lAikawall2007l) . 
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On the other hand, lOberg et aTl (l2007l) measured the pho- 
todesorption rate of CO in the laboratory. They found a rate sim- 
ilar to that found for H2O by Westlev et al. ( 1995), about 2 orders 
of magn itude higher than prev iously available theoretical esti- 
mates ( Draine & SalpeteJl 19791) . implying that this non-thermal 
desorption process was an appealing candidate for explaining the 
abundance of cold CO in various astrophysical environments. 
This motivated us to investigate the influence of both vertical 
mixing and photodesorption on the CO column density in disks, 
using a parametric model consistent with interferometric obser- 
vations of circumstellar disks. 

In Sect. 12] we describe the physical and chemical model 
used, while Sect.|3]is devoted to the results of our computations. 
In Sects. |4]and|5] we discuss our results and present our conclu- 



2. Modeling 

2.1. Disk structure 

To consistently co mpare with the physical structure observed by 
jPiefii et all ( l2007h . we use a simple disk model with physical 
quantities varying as power la ws of the radius, similar to the 
parametric disk model used by iPietu et al. (2007) to reproduce 
their interferometric maps. The vertical temperature distiibution 
is defined to have two layers, one cold layer with a constant tem- 
perature of 10 K, and a warm layer in radiative equilibrium with 
the central star radiation. This behavior is typical of disk mod- 
els (e.g. D'Alessio et al. 1999; Aikawa & Nomura 2006). The 
precise location of the transition between the two layers varies 
between models. We chose this transition at 2 pressure scale 
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heights, a location similar to that used by iDartois et alj (l2003h 
to fit their data. The constant mid-plane temperature, although 
surprising at first sight, is consistent bot h with the temperature 
observed in '^CO J=l-0 bv iPietu et all (l2007i) and the surface 
density power exponent of -3 /2 observed in the outer parts of 
most "un disturbed" disks (see DutrevllToOSl for a review). In an 
a-model ("Shakura & S vunvaevll973l K the product of the temper- 
ature and the surface density is indeed assumed to vary as R"^^^^. 
Therefore, a constant mid-plane temperature in an stationary a- 
model is the only way that the surface density decreases as R"^^^. 
This regions of almost constant temperatu re are also found in the 
more complete disk models developed bv iP' Alessio et al.] (1 19991 
I200I . 



The laws that we use in the present paper are: 
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where Cj is the isothermal sound speed, Q is the Keplerian fre- 
quency, and E is the surface density. From these laws, the mass 
density p is then computed from the hydrostatic equilibrium. 

2.2. Chemistry 

The chemical evolution of the disk is computed with a ID verti- 
cal model that takes into account both gas-phase and grain sur- 
face chemistries. This code was adapted from the OSU gas-grain 
code developed by the astrochemistry team of Eric Herbst. The 
formalism of surface chemistry and various processes a re de- 
scribed in lHasegawa et all (1 1 9921) and lGarrod et al.l (12007 ). The 
gas-phase chemical network uses updates from osu_83_2O0^. 
The complete network (gas and grain) contains 655 species and 
6067 reactions. Grains are assumed to be of interstellar type. To 
model the C O and H2 UV self-shielding, we use the approx- 
imation from iLee et al .' (1996). The code is written in Fortran 
90 a nd uses the DLSO DES routine from the ODEPACK pack- 
age dHindmarshl 1 1 9831) to solve the differential equations. The 
initial abundances of our calculations are taken from the out- 
put of a gas-phase calculation of a dark cloud (nH=2 lO'* cm"-', 
T - 10 K, Av= 10) after 1 0^ yr and our eleme ntal abundances 
are the EA3 conditions from lWakelam & Herbst. (2008), with an 
S abundance of 8 x 10 The gas-grain model is then allowed 
to evolve for 10^ yr with 32 vertical points. A period of 10^ yr 
is the radial mixing timescale at 100 AU, as well as the evolu- 
tionary timescale for the disk structure. Computing for a longer 
time would in principle require a full 2D (R,z) approach and an 
evolving disk model. 



2.3. Photodesorption 

In their experiments, 'Ober g et"an (l2007h found a photodesorp- 
tion rate Y po of 3 x 10^^ CO molecules per UV photon above 
pure CO ice. In circumstellar disks, UV photons originate from 
two main processes: (i) the ambient UV field, coming from both 



the central star and the surrounding ISM, and (ii) secondary pho- 
tons generated by the interaction of H2 with cosmic rays and cos- 
mic ray-induced (CR) secondary electrons (Prasad & Tarafd;^ 
ll983l:IShen et al.l2004j) . Photodesorption of CO (we did not con- 
sider photodesorption for other species) by these two sources of 
UV photons can be expressed by the chemical reactions: 



adsorbed CO 
adsorbed CO 
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where kuv a n d kcu are the first-order rate coefficients (s ' ). From 
lOberg et al.l (I2OO7I) and the additional approximation that ad- 
sorbed CO molecules can always escape independently of the 
thickness of the adsorbed layer, we can express these rate coef- 
ficients as: 
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(Hassel et al. in preparation) where Iism- 
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10** photons 

cm"^ s ' is the a mbient far UV field strength for a standard ISM 
irradiation field (lDrainelll978h . 

IcR-Fuv - 10^ photons cm ^ 
s ' is the strength of the far UV field generated by the cosmic 
rays, y measures the difference between the UV extinction and 
the vi sual extinction {y ~ 2 for interstellar grains. iRoberge et"an 
1991), and cTsiies is the site density per grain surface unit (i.e. 
the maximum number of molecules adsorbed as a single layer 
on a grain per cm^). The quantity fuv is a scale factor used to 
express the strength of the ambient UV field (originating mainly 
in the central star) as a multiple of the standard ISM radiation 
field. This expression is justified by noting that the UV spectrum 
coming from the UV excess of young stars has a similar sh ape 
to that of the standard interstellar field dChapillon et al.l2008h . In 
the present report, we use: 
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The Uy field intensity is consistent with the results of 
iBergin et alj (12003). We consider here only the vertical extinc- 
tion of the UV field. This approximation is motivated by the fact 
that the direct (radial) extinction of stellar UV radiation is very 
efficient so that most photons seen by the disk are in fact scat- 
tered by small grains at high altitudes. We consider here that this 
scattering acts at the edge of our computation box (4 pressure 
scale heights). With this law for fuv, direct photodesorption is 
dominant for Av < 7, while photodesorption by CR-induced UV 
photons (independent of Av) maintains a low desorption rate in 
denser regions. 

2.4. Turbulent mixing 

We comp ute the vertical turbulent mi xing using an a turbulent 
viscosity dShakura & Svunvaevlll973t) . where the diffiisivity D, 
is given by: 



Q- Sc — 
Q 



(11) 



The Schmidt number Sc, defined as the ratio of the turbulent 
viscosity to the turbulent diffiisivity, is a poorly constrained 
parameter. In the present work, we set its val ue to be 1 and 
choose a - 10"^, values c o mmonly used ( e.g. ISemenov et al] 
(l2006l) : IWillacv et all 120061) : lAikawal(l2007h '). 
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Fig. 1. Vertical distribution of gaseous CO at 100 AU from the 
central star, (a) no photodesorption, no diffusion (dotted-dashed 
line), (b) no photodesorption, diffusion (dashed line), (c) pho- 
todesorption, no diffusion (dotted line), (d) photodesorption, dif- 
fusion (solid line). Superimposed on the graph are the tempera- 
ture (thick grey line) and visual extinction (thick dashed grey 
line) as a function of height above the midplane. 

The resulting equations for the chemical abundances x, are: 

^^P,-Li + --D,p-Xi (12) 
at p oz oz 

where f , and L, are the chemical production and loss terms. 
The vertical mixing and the chemistry are computed using first- 
order operator splitting, the chemistry being computed before 
the mixing inside a time step. This induces intrinsic errors, which 
are controlled by reducing the time step. We adapted the time 
step so that the splitting errors remain below 5 % relative accu- 
racy, by comparing with simulations without operator splitting, 
which are far more time-consuming. 

3. Results 

The vertical distribution of gas phase CO was computed at 100, 
200, 300, 400, 500, and 600 AU, for 4 different models, with 
or without CO photodesorption and/or vertical mixing. Figures 
[1] |2] and [3] show the resulting vertical CO distributions. The 
"handle-bar" shape of the vertical abundance profile at all radii 
is the result of CO photodissociation at high altitudes (mostly 
visible at 300 and 500 AU) and CO freezing onto grains in the 
dense mid-plane (below z/// - 2). Carbon is locked in different 
dominant forms as a function of height; frozen CO in the mid- 
plane, then frozen CO2, gaseous CO, and finally atomic C in the 
photodissociated region. The thickness of each layer depends on 
the radius (through density, temperature, and Av). In the 100 AU 
profile (Fig. [TJ, a dip is visible in the CO distribution without 
mixing. This dip is the consequence of a layer of N2H^, which 
initiates a secondary reaction chain (dominated by grain chem- 
istry) started by: 

N2H+ + CO ^ N2 + HCO+. (13) 

HCO'*^ is then converted into other carbon-bearing 
molecules, the reaction chain ending with frozen CO2. 
The CO column density as a function of radius is displayed in 

Fig. a 
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Fig. 2. As Fig. 1 but for R=300 AU 
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Fig. 3. As Fig. 1 but for R=500 AU 



3.1. Influence of photodesorption 

Photodesorption by primary UV photons is the most important 
process affecting the abundance of gas-phase CO. When the vi- 
sual extinction is sufficiently low (below about 5 mag), this des- 
orption mechanism is efficient enough to preserve a significant 
abundance of CO. For radii larger than 200 AU, the total column 
density of gas phase CO is increased by an order of magnitude 
when this process is taken into account. The situation is differ- 
ent at 100 AU. Here, the visual extinction reaches a value of 5 
mag at 2 pressure scale heights, and grains receive insufficient 
primary UV photons to enable efficient photodesorption of CO. 
The secondary UV photons originating from the interaction of 
H2 with electrons induced by cosmic rays are then the only pho- 
todesorbing photons. They maintain a small fraction of CO in 
the gas phase but do not influence the total CO column density 
by more than 10 %. 

3.2. Influence of vertical mixing 

In our simulations, vertical mixing has a small influence on the 
total column densities of CO. At 100 AU (and 200 AU), as 
shown in Fig. H) the presence of a layer of N2H^ appears as a 
sink of CO (via reaction [13) and the total CO column density is 
reduced by vertical mixing. This is the reason why the CO abun- 
dance below z - 2H is lower in the case with diffusion (Fig.lTJ. 



4 



F. Hersant, V. Wakelam, A. Dutrey, S. Guilloteau and E. Herbst: Cold CO in circumstellar disks 



E 











- 










_ _ — 






/ 






<' ; 






/' : 






// 






// 






//' 






//■ 1 





















o 



nOO 200 500 1000 

Radius(AU) 

Fig. 4. Total column density of gaseous CO as a function of 
the radial distance, (a) no photodesorption, no diffusion (dotted 
dashed line), (b) no photodesorption, diffusion (dashed line), (c) 
photodesorption, no diffusion (dotted line), (d) photodesorption, 
diffusion (solid line). The grey lines are the corresponding cold 
CO (T < 17 K) column density profiles. 

The low efficiency of vertical mixing is a consequence of two 
effects. First, when CO condenses onto grains, grain chemistry 
transforms it into more refractory species, thereby preventing 
desorption when it diffuses back to the warm region. The sec- 
ond reason is the consequence of a dilution effect. Due to a tem- 
perature transition at 2 pressure scale heights, the warm region 
represents less than 10 % of the total mass. This is obviously a 
possible bias of our disk model, the influence of vertical mixing 
being very sensitive to this transi tion. For a transition deepe r in- 
side in t h e disk , as considered bv lAikawa & Nomural (l2006l) and 
lAikawal (|2007|) . vertical mixing is far more efficient. Accurate 
knowledge of the true thermal structure of the disk would be re- 
quired to obtain a more precise estimate. 

4. Discussion 

Figure ID shows that photodesorption increases the total CO col- 
umn densi t y, but also makes the CO column density flatter 
iPietu et al] (|2007|) measured a steep radial profile of CO column 
density, varying as R"^ or even more steeply, with about 10'^ 
cm"^ at 300 AU. This gradient is steeper than for our models 
without photodesorption. Quantitative comparisons with obser- 
vations are however far from being straightforward because they 
require the generation of CO emission maps and iterations on 
the physical structure to find the best fit for each chemical setup. 
Incidentally, the flat CO column density profile is a common fea- 
ture of disk chemical models. In our case, we would like to raise 
a few points: 

- Our results are extremely sensitive to the location of the tran- 

sition between the cold and warm regions. We chose a tran- 
sition at 2 pressure scale heights; other choices would lead to 
significant differences in the column density radial profile. 

- At large distances from the central star, the CO column density 

is very sensitive to photodissociation. The stellar UV flux 
is not constrained accurately, although its value can signifi- 
cantly alter our results. Our treatment of the UV field is also 
approximate. We only consider UV radiation in a direction 
perpendicular to the disk. Observations of ' ^CO are sensitive 
only to the outer regions jPietu et al.ll2007l) . where the inter- 



stellar UV field impinging from the disk edge will steepen 
the slope of the CO column density profile. 

- For a given UV flux, photodissociation is controlled by the re- 

lation between visual extinction and column density. This 
factor depends mostly upon the amount of small grains. 
Grain growth can change the penetration of UV radiation in- 
side the disk and dissociate CO further (see IChapillon et al.l 
|2008). Since grain growth rates depend upon density, the 
grain distribution can depend upon the radial distance and 
completely change the radial profile of the CO column den- 
sity. In turn, it may be possible to invert the CO column den- 
sity to infer some information about the current state of grain 
growth at a given radius. 

- Grain growth also increases the timescales for the adsorption 

of CO and grain surface chemical reactions, making both 
photodesorption and vertical mixing more efficient at retain- 
ing CO in the gas phase. 

5. Conclusions 

We have shown that CO photodesorption appears to be a good 
candidate to explain the abundance of cold gas phase CO in cir- 
cumstellar disks. However, it tends to create a flatter radial pro- 
file for the CO column density, in contrast to that observed in 
disks. Therefore, more work is necessary to solve this problem. 
Consideration of the photodissociating UV flux and variations 
in its penetration inside the disk due to grain growth and geom- 
etry seem to provide an avenue to improve our modelling of cir- 
cumstellar disks. Moreover, we have shown that the influence of 
vertical mixing on the abundance of CO is not straightforward. 
Rather, the efficiency of vertical mixing depends upon many pa- 
rameters, including the details of the disk structure (density and 
temperature) and the level of grain growth, all of which remain 
uncertain today. 
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